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Abstract 

We consider an integrable hierarchy of nonlinear evolution equations (NLEE) re- 
lated to linear bundle Lax operator L. The Lax representation is 7Li x Z2 reduced 
and is naturally associated with the symmetric space SU(3)/S(U(1) x U{2)). The 
simplest nontrivial equation in the hierarchy is a generalization of Heisenberg ferro- 
magnetic model. We construct the iV-soliton solutions for an arbitrary member of 
the hierarchy by using the Zakharov-Shabat dressing method with an appropriately 
chosen dressing factor. Two types of soliton solutions: quadruplet and doublet soli- 
tons are found. The one-soliton solutions of NLEEs with even and odd dispersion 
laws have different properties. In particular, the one-soliton solutions for NLEEs 
with even dispersion laws are not traveling waves; their velocities and their ampli- 
tudes are time dependent. Calculating the asymptotics of the iV-soliton solutions 
for t — > ±00 we analyze the interactions of quadruplet solitons. 
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1 Introduction 

5— 1 

The main object of present paper is the following coupled system of equations 

iu t + u xx + (uu* x + vv*)u x + {uu* x + vv*) x u = 
iv t + v xx + [uu* x + vv*)v x + [uu* x + vv* x ) x v = ' 



where the smooth functions u : M 2 — > C and v : IR 2 — > C satisfy the algebraic constraint 
\u\ 2 + \v\ 2 = 1. The system ( 11 .11) is a natural candidate to be a multicomponent general- 
isation of the classical Heisenberg ferromagnetic equation. It is well known [32J that the 
Heisenberg ferromagnetic model is integrable in the sense of inverse scattering method 
(ISM). It has a Lax pair related to the algebra su(2). Since the time the complete inte- 
grability of HF equations was discovered, many attempts for its generalization have been 
made |2Ql EH E2j ■ A well known method [26l EZl EHl [12l M IM1 EOl [31] to obtain new 
integrable nonlinear evolution equations (NLEE) is based on imposing certain algebraic 
reductions on generic Lax operators. Lax pairs associated to hermitian symmetric spaces 
represent a special interest in modern theory of integrable systems is study of NLEEs 
PO El EH] since the NLEEs they produce look relatively simple. 
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The system (II. ip is also integrable in the sense of ISM. Its Lax operators are associated 
with the symmetric space SU(3)/S(U(1) x U{2)) with a Z 2 x Z 2 reduction imposed on 
them [I5l [161 H3]. 

The purpose of the present paper is to derive the soliton solutions for the integrable 
hierarchy of equations related to (11. ip and analyse the interactions between them. In this 
sense this paper is a natural continuation of our previous papers [15j [161 fl3] . 

In Section 2 we start with some basic facts to be used further in paper. Firstly 
we describe the hierarchy of nonlinear equations related to (11. ip in terms of recursion 
operators. Then we outline the spectral properties of the relevant Lax operator and 
formulate direct scattering problem. The spectrum of scattering operator L consists of 
a continuous and a discrete parts. As a result of the Z 2 reductions L possesses two 
configurations of discrete eigenvalues: generic ones, coming in quadruplets ±A/c, ±A£ and 
purely imaginary ones coming as doublets ±ittj. 

In Section 3 we derive the 1-soliton solutions for the NLEEs of the hierarchy. For 
this to be done we apply the Zakharov-Shabat dressing method [34J [38j |35j [36] with a 
rational dressing factor with 2 simple poles. Due to the action of reductions we have 
two types of 1-soliton solutions: quadruplet solitons to correspond to 4 eigenvalues and 
doublet ones to correspond to 2 eigenvalues respectively. We present explicit expressions 
for these two types of one-soliton solutions. In order to construct general multisoliton 
solutions we discuss two different purely algebraic constructions: by using a multiple pole 
dressing factor and by applying "one-soliton" dressing factors several times consecutively. 
It turns out that the properties of the 1-soliton solutions to NLEEs with even and odd 
dispersion laws differ drastically. For example, the 1-soliton solutions for NLEEs with even 
dispersion laws are not traveling waves. Even the doublet soliton of eq. (II .ip exhibits two 
maxima (resp. minima) for (resp. for \vi\) which first come closer to each other and 
then move away, one to oo and the other to — oo as time goes to t — > oo. Their velocity, 
as well as their amplitudes are time dependent. These properties are similar to the ones 
of the boomerons and trappons discovered by Calogero and Degasperis j2[ [31 HJ [5] . At the 
same time the soliton solutions to the NLEEs with odd dispersion laws (e.g. the solutions 
of eq. (I2.18P ) behave as standard solitons, i.e. they are traveling waves. 

Section 4 is dedicated to interactions of quadruplet soliton solutions for the NLEE 
with odd dispersion laws. In order to do this we use the classical method of Zakharov 
and Shabat, see the monographs [34], |32] for a detailed exposition. Namely, we calculate 
the limits of the iV-soliton solutions for t — > ±oo assuming that all solitons move with 
different velocities. In this way we establish that the solitons preserve their velocities and 
amplitudes; the only effect of their interaction consists in shifts of the relative mass center 
and the phase of solitons. We provide explicit expressions for these shifts in terms of the 
poles Hk of the dressing factors. 

In Section 5 we briefly discuss the conservation laws of the NLEE and finish with some 
conclusions. 

2 Preliminaries 

In this section we shall expose in brief some basic facts on Lax operators and direct 
scattering problem for the integrable hierarchy of the equation (II. ip . In doing this we 
shall use a gauge covariant formulation [TH [TTJ, HH [19]. 
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2.1 Polynomial Lax Pair Related to SU(3)/S(U(1) X U(2)) 

The NLEEs under consideration in this paper represent a zero curvature condition [I, A] = 
for Lax operators L and A in the form: 

L(X) = id x + XLi(x,t) (2.1) 

N 

A(A) = i^ + ^A fc A fc (x,t), (2.2) 

k=l 

where A G C is the so-called spectral parameter and the functions I\ and A^, k — 1, . . . , N 
take values in sl(3, C). The Lax operators are subject to the following Z 2 reductions: 

L f (A*) = —L(X), A\\*) = -A{\) (2.3) 
CL(-A)C = L(A), CA(—X)C = A(X), (2.4) 

where C = diag (1,-1,-1) and the operation "is defined as follows 

L(X)ip(x, t, A) = id x i[)(x, t, A) — Xi[)(x, t, X)Li(x, t, A). 

Due to reduction (12. 3p the matrix coefficients of the Lax pair are hermitian matrices. On 
the other hand reduction (12 .4p represents an action of Cartan's involutive automorphism 
which defines the symmetric space SU(3)/S(U(1) x U(2)), see [221 It induces a 

Z 2 -grading in the Lie algebra sl(3, C) 

sl(3) =sl°(3) ®8l\3), sf (3) = {X gs((3) | CXC = {-l) a X}. (2.5) 

It is evident that Li,Ak G si (3) for being an odd integer and G sl°(3) otherwise. 
This means that A^ for even are block-diagonal matrices of the form 



-4 A 



while Li and A^ for odd k have the complementary block structure. In particular, L\ is 
written as: 

/ u v\ 

L x = I u* (2.6) 
\v* 0/ 

The potential Li is required to obey the following conditions: 

1. The eigenvalues of L\ are 0, ±1, i.e. the potential satisfies the characteristic equation 
L\ = In- 




2. The function Li(x,t) — L± where 



e^± 



lim L x (x,t) = L±= \ 0, <f>± G R. (2.7) 

1 e - 

is a Schwartz type function, i.e. it is infinitely smooth and tends to faster than 
any polynomial when \x\ — > 00. 
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The grading ( 12. 5 p means that any function X with values in sl(3) can be split as 
follows: 

X = X° + X\ X ' 1 G s( 0,1 (3). (2.8) 

Let us define the Killing form for s((3) as follows: 

(X,Y) =tr(XY), X,FGsl(3). 

Then each component X ' 1 splits into a term commuting with L\ and its orthogonal 
complement with respect to the Killing form 

X° = X ^ + k L 2j L 2 = L\- jjl, (X ^, L 2 ) = (2.9) 
X 1 = X 1 ^ + kjLi, (X 1 ' ± ,L 1 ) = 0. (2.10) 
As a simple consequence of condition 1 above L\ and L 2 are normalized as follows: 

(L U L 1 ) = 2, (L 2 ,L 2 > = | (2.11) 

Therefore the coefficients and K\ are given by the following equalities 

k = ^(X°,L 2 ), Kx = ±(X\L x ). (2.12) 

The zero curvature condition [L,A] = for the pair (12. ip . (12. 2p leads to certain 
recurrence relations for the matrix coefficients of L and A, see [13]. Resolving them 
allows one to express Ak in terms of L\ and its x- derivatives of order up to X — k. Since 
the maximal order term in the operator A must commute with L\ there exists two options: 

a) A N = c 2p L 2} if N = 2p 

b) A N = c 2p+1 L u if N = 2p + 1, 

where c 2p and c 2p+x are constants. Then a more detailed analysis [13] shows that the 
NLEEs look as follows: 

v p-i 

a) ad^Li,* - ^ c^AiA^^ad^ 1 ^ - ^ c 2g+ i (AiA 2 ) 9 ad l\L ltX = 0, 

9=1 g=0 
V V 

b) ad^Li.t - ^c 2f ,(A 1 A 2 ) 9 - 1 A 1 ad^ i 1 L 2 , ;c - ^ c^AiA^ad^ 1 !^ = 0. 

9=1 9=0 

The integro-differential operators Ai and A 2 appeared above are given by 

Ai = -iad^ 1 (W) - ^d'^d^L,)) 

A 2 = -iad^ 1 (ird x (.)- h^d' 1 (d x (-), L 2 ) 

where projection n := ad^ad ^ cuts all Lx-commuting parts off. The operator 

(A X A 2 X, Xgs1°(3) 
\ A 2 AxX, X gs[ x (3) 

is called recursion operator. It can be viewed as an adjoint representation of the operator 
L. Its existence manifests the hierarchies associated with NLEE (nonlinear equations, 
integrals of motion, simplectic forms etc) and thus plays a very important role in theory 
of solitons. 



(2.13) 



(2.14) 
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Example 1 Consider the simplest case when N = 2. Then the matrix coefficients of the 
second Lax operator A read: 

/ 1/3 \ / a b\ 

A 2 = - \u\ 2 -2/3 u*v , Ai=(a*00] (2.15) 
\ v*u \v\ 2 -2/3 ) \b* 0/ 

a = \u x + \{uu* x + v v l)u, b = iv x + i{uu* x + v vl)v (2-16) 
This L-A pair produces the 2-component system 



iu t + u xx + (uu* x + vv* x )u x + (mm* + vv* x ) x u = 

IVt + V xx + (MM* + W*)Ua: + + W*)^ = ' 



(2.17) 



we 



started our paper with (see ( fi.ij) ).D 



For completeness here we present another member of the hierarchy (12. 13ft . It is the 
simplest NLEE corresponding to an odd dispersion law. 

Example 2 Consider the case when f(A) = — 8A 3 J ; i.e. C3 = —8, C2 = c\ = 0. Then the 
corresponding 2-component system obtains the form: 



u t = Su xxx + 12(uu* x + vv*)u xx + r(u, v)u x + s(u, v)u 
v t = 8v xxx + 12(uu* x + vv* x )v xx + r(u, v)v x + s(u, v)v ' 



(2.18) 



where 



r(u, v) = 3 [4:(\u x \ 2 + \v x \ 2 ) + 5(uu* x + vv* x ) 2 + Q(uu* x + vv* x ) x ] 

s(u, v) = 3 [2(uu* x + vv* x ) xx + 4(\u x \ 2 + \v x \ 2 ) x + 5(uu* x + vv* x ) 2 x ] . □ 

Sometimes it is more convenient to deal with Lax operators written in canonical gauge. 
In this gauge the operator (12. ip looks as follows: 

L(\) = g- 1 Lg = id x + U (x,t) + \J, J = diag (1, 0, -1), (2.19) 

where 

f2 ( 1 ° _1 ' 

g = y-\u* V2v u* ] . (2.20) 

\v* -V2u v* 
The second Lax operator (12.21) is given by 

AT-l 



a) 1(A) = id t + J2 xk Mx,t) + c N X N I, N = 2p 

k=0 
N-1 

b) i(A) = id t + ^\ k A k {x,t) + c N X N J, N = 2p+1, 
where I = g~ l L 2 g = diag (1/3, -2/3, 1/3). 



(2.21) 
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2.2 Direct Scattering Problem 

In order to formulate a direct scattering problem for L, one needs to introduce auxiliary 
spectral linear system 

L(X)if)(x, t, A) = id x if)(x, t, A) + ALiO, t)if)(x, t, A) = 0. (2.22) 

Here if) denotes a fundamental set of solutions or a fundamental solution for short. Since 
the operators (I2.ip and (12.21) commute if> also satisfies 

A(X)if;(x, t, A) = \ idt + xkA k( x , *) J ^fo t, A) = ip(x, t, A)f (A) (2.23) 
as well. The matrix- valued function 

N 

f(A)= lim fo'VA^fMlfe (2.24) 

x— >±oo ' 4 

h=l 

is called dispersion law of the nonlinear equation (12.131) . The unitary matrix 



g ± = lim g(x,t) = -±= | V2e 1 ^ 

X — >±oo 2, 



1 -1 

D s/2^± 



involved in the definition of the dispersion law diagonalizes the asymptotics L^± = 
lim^-too Li(x, t). It can be proven that the dispersion law of (I2.13P reads 

p-i p 

a) f (A) = C 2q+l^ 2q+1 J + C 2^ 29/ ' 

9=0 9=1 (2.25) 

b) f (A) = J2 C2 q+ i\ 2q+l j + Yl c ^ x29L 

q=0 q=l 

The dispersion law of the 2-component system (12.171) is —A 2 / and that of ( (I2.18P is — 8A 3 J. 
It is evident from (I2.25P that f(A) obeys the splitting: 

f(A) = f (A)/ + f 1 (A)J, (2.26) 

which is a result of the Z 2 grading (12.51) of the algebra sl(3). 

A special type of fundamental solutions are the so-called Jost solutions if>± which are 
normalized as follows 

lim ijj ± (x,t,X)e- iXJx g ± 1 = 1. (2.27) 

x— >±oo 

Due to (12.241) one can show that the asymptotic behavior of if>± do not depend on time 
and thus the definition is correct. The transition matrix 

T(t, A) = [^+(x,t,A)]-V-(^t,A). (2.28) 

is called scattering matrix. It can be easily deduced from relation (I2.23P that the scattering 
matrix evolves with time according to the linear differential equation 

id t T + [f(X),T) =0, (2.29) 
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which is integrated straight away to give 

T(t, A) = e if(A) *T(0, A)e" if(A) *. (2.30) 

From now on the parameter t will be fixed and we shall omit it to simplify our notation. 
Due to reasons of simplicity we set <p + = 0_ = as well. 

The action of Z2-reductions ( 12. 3ft . ( 12.4ft imposes the following restrictions 

V£(*, A*)] 1 = if, ± (x, A), [Tt(A*)] - 1 = T(A) 
C^±(x,-A)C =V±(ac,A), CT(-A)C=T(A) 

on the Jost solutions and the scattering matrix. 

The continuous spectrum of L fills up the real axis in the complex A-plane. Thus the 
A-plane is divided into two regions denoted by C + (the upper half plane) and C_ (the 
lower half plane). These regions represent domains for fundamental solutions x + {x,X) 
and x~( x ,X) to be analytic functions in C + and C_ respectively |16j . The fundamental 
analytic solutions (FAS) can be constructed by using Gauss factors in the decomposition 
of the scattering matrix: 

T(A) = T T (X)D ± (X)(S ± (X))- 1 . (2.32) 

S + and T + are upper triangular matrices, S~ and T~ are lower triangular matrices and 
are diagonal ones. Then x + an d X~ are expressed as follows 

X ± (x, A) = ^(x, X)S ± (X) = ^ + (x, X)T^(X)D ± (X). (2.33) 

Due to relation (12. 33ft the FAS can be interpreted as solutions to a local Riemann-Hilbert 
problem 

X + (x, A) = X -(x, X)G(x, A), G(X) = (S~ (X))- 1 S+ (X) . (2.34) 

The established interrelation between the inverse scattering method and Riemann-Hilbert 
problem plays an important role in constructing solutions to NLEEs through dressing 
method. 

It can be shown that the reduction conditions (12.311) and equation (12.321) lead to the 
following demands on the Gauss factors 

[S + (A*)] f = [S-(A)]-\ CS ± (-X)C = S T (X) 

[T+(A*)] t = [T-(A)]"\ CT^-AJC = T T (A) (2.35) 

[D+(X*)] ] = [D-(X)]-\ CD±(-X)C = D±{\), 



where 



C 




Finally, combining all this information we see that the FAS obey the symmetry conditions 
[ X + (x, A*)] = [ X -(x, A)]" 1 C X + (x, -A)C = X ~(x, A). (2.36) 
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3 Dressing Method and Soliton Solutions 



As we mentioned in the previous section the inverse scattering method is tightly related to 
Riemann-Hilbert problem. The Riemann-Hilbert problem possesses two types of solutions: 
regular ones (without singularities) and singular ones. Singular solutions can be generated 
by dressing regular solutions with a factor which has prescribed singularities. The simplest 
types of singularities are first order poles and zeroes. It can be proven that they correspond 
to poles of the resolvent of L. Hence they are discrete eigenvalues of the Lax operator 
(12. ip . The discrete eigenvalues of L form orbits of the reduction group Z 2 x Z 2 . There 
exist two types of orbits: generic orbits containing quadruplets of eigenvalues {±/x, ±//*} 
and degenerate orbits consisting of two imaginary eigenvalues ±i/t (doublets). 

There is a very deep connection between singular solutions to Riemann-Hilbert prob- 
lem and soliton solutions to the corresponding nonlinear problem. In the present section 
we are going to analyze the soliton solutions to the system (I2.13p . For this to be done, we 
are going to apply the dressing method proposed in [38] and developed in (351 EHl l28| 129]. 
We demonstrate that the NLEE (I2.13P has two types of 1-soliton solutions: doublet soli- 
ton to be connected with two imaginary discrete eigenvalues of L and quadruplet soliton 
connected to 4 eigenvalues. 



3.1 Rational Dressing 

The dressing method is an indirect method for solving a NLEE possessing a Lax represen- 
tation. This means that it allows one to generate a solution to the NLEE starting from 
a known one. Let us assume we know a solution 

/ u v ' 
4 0) = \ u* 
\v* 

of (12.131) and a fundamental solution ipo(x,t, A) of the auxiliary linear problems 

L(°)(A)^ = i^o + ALfVo = 

N (3 1) 

A(°)(A)^o = ifyfo + Yl = 0. 

fc=i 

Then one constructs another function ipi(x, t, A) = $(x, t, \)ipo(x, t, A)$1(A), where $-(A) 
lim^-oo $(x, t, A). This function is a common solution to 

L«(A)Vi =i9,^i + AL[ 1 Vi =0 

N (3 2) 

a«(a)^ = ic^i + = o, 

fc=l 



where the potential 



U\ V\ 
L[ 1] = ( u\ 
v* x 



s 



is to be found. From (13. ip and (13. 2p it follows that the dressing factor $(a;, t, A) satisfies 
the following equations: 

id x $ + XL^h - A$4 0) = (3.3) 

TV N 

idt$ + '%2\ k A < £ ) $-$'22\ k A ( £ ) = 0. (3.4) 

fc=l k=l 

We also assume that the dressing factor is regular at |A| — > 0, oo. Then from (I3.3P one 
can derive the following relation between Lp and L^: 

L { i\x,t) = $(x,t,oo)Lf\x,t)&(x,t,oo). (3.5) 

This equation will play a central role in our further considerations since it allows one to 
generate a new solution to ( 12. 131) from the given one . 

Due to the reduction conditions (1 2 . 3 j) . (12.41) the dressing factor obeys the symmetries: 

C<f>{x,t,-X)C = ${x,t,X) (3.6) 
$(x,t,\)&(x,t,\*) = H. (3.7) 

In order to obtain a nontrivial dressing we choose <&(x, t, A) as a rational function^ of 
A with a minimal number of simple poles. At first we shall consider the case when these 
poles are generic complex numbers. Hence the dressing factor looks as follows: 

n , XM(x,t) XCM(x,t)C /ooN 

${x, t, A) = 11 + 1 ; + ' , (3.8) 

A — (i A + [i 

where Re \i ^ 0, Im fi ^ 0. It is evident that the reduction condition (I3.6P is fulfilled. On 
the other hand (13. 7|) leads to the conclusion that 

t -(x,u) = a + Wtfet) + AC f (l ' t)C - (3-9) 

A — n* A + /i* 

The identity $(A)$ _1 (A) = 11 must hold for any A. Therefore after equating the 
residue at A = fi* to one gets the equation: 

„ u*M(x,t) u*CM(x,t)C\ ,, t/ . 

i + a — k —lJ- + n v ' ; mUx, t) = o. 3.10 

A** — A* /i* + A* / 

The rest of algebraic relations can be reduced to (13.101) due to the symmetry conditions 

(Q, (1231). 

The residue M ought to be singular since otherwise it should be proportional to 11 and 
the dressing becomes trivial. It suffices to consider the case rankM = 1. Then M can be 
decomposed in the following manner: 

M=\n)(m\, \n) = (m, ri2, ns) T , (m\ = (m*, m*,, rn,g). (3-H) 

After substituting this representation into (13. lQf) one derives a linear system for the 3- 
vector \n): 

H*\n){m\m) fi*C\n)(m\C\m) _ 
lm) 2i^ + %j " U - [6AZ) 



1 If $ is A-independent then it does not depend on x and t either. Thus (|3.5[) produces simply a unitary 
transformation of Lj° which is not essential because of U(2) gauge symmetry of the model. 
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where we have used the notation u — Re (J,, k — Im (i. The solution of (I3.12p reads: 



, . 1 f (m\m) (m\C\m) \ . , , 
\n) = — i-M - ' ' C m). 3.13 
1 ' /i* V 2ik 2w / 

The vector |m) is an element of the projective space CP 2 . Indeed, it is evident that a 
rescaling |m) — > h\m) with any complex h ^ does not change the matrix M. 
Taking into account the ansatz (13.81) one can rewrite (I3.5P as: 

4 1} = (1 + M+ CMC)4 0) (11 + M + CMC) 1 . (3.14) 

Notice that the dressing procedure preserves the matrix structure of L since the factor 
11 + M + CMC is a block-diagonal matrix. 

We have expressed all quantities needed in terms of \m) and now it remains to find 
\m) itself. For that purpose we rewrite equations (13.31) . (13.41) in the form: 

$(s, t, a) (ia, + al! 0) ) q>-\x, t, a) = alS 1} 

/ N \ N (3.15) 

9( X , t, a) [id t +j2 A fc 4 0) < ^ 1 (^> t,X) = J2 A fe 4 1} - 

V k=l J k=l 

It is obviously satisfied at A = 0. After equating the residues of (13.151) at A = ji* to we 
obtain a set of differential equations 

/ /i*M /i*CMC ^ / (0) n 
II + — + — — I id x + /i L\ ) \m) = 

Taking into account (13.101) the equations above can be reduced to 

id x + (i*L^ (x, t)j \m(x, t)) = h(x, t) \m(x, t)) 

(3.17) 



(3.16) 



N 



id t + J2^*) kA k\ x ^) \m(x,t)) = h(x,t)\m(x,t)) 



k=l 



for some arbitrary function h. At this point we recall that the vectors in the decomposition 
(I3.1ip are not uniquely determined. Indeed, the operation |n) — > B~ l \n) and |m) — > B*\m) 
for any nondegenerate 3x3 matrix B produces another decomposition of M. It is not 
hard to see that it is always possible to choose B in such a way that h = is fulfilled. 
Thus from (I3.17P it follows that \m(x,t)) is proportional to some fundamental solution 
ipo(x,t, A) of the bare linear problem, namely 

\m(x,t)) = tp (x,t,[i*)\m ), (3.18) 

where |m ) G C 3 \{0} is a constant vector of integration. The new solution L{ 1] of ([235} 
and the solution ^(x, t, A) of the corresponding linear system are parameterized by a 
complex number (i and a complex 3- vector \mo)- 
Thus we have proved the following Proposition: 
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Proposition 1 Let Lj ' be a solution of A2.13\) and ijjo(x,t,X) be a common solution to 
Ii3. Let also \i be a complex number to fulfill Re \x ^ 0, Im /i > and |m ) G C 3 \{0}. 
Then the matrix-valued function L^\x,t) defined by ([3.141) where M = \n)(m\ is deter- 
mined by A3.13A) and A3. ISA) is a solution to A2.13\) as well. The corresponding fundamental 
solution ipxfajt, A) of A3. 2]) is given by ipi = Qipo where ®(x,t, A) is determined by A3. S]) . 
AJU\) . ATT3}) and AJ7TS]) . □ 

Let us now consider the case when the poles of the dressing factor are imaginary, i.e. 
we have: 

r M(x,t) CM(x,t)C 



<S>(x,t,\) = Z+\(-^+ K [ J , k^O. (3.19) 

\ A — IK A + IK J 

Then <3> _1 has the same poles as $ and therefore the equality = 11 already con- 

tains second order poles. In this case the natural requirement of vanishing of the matrix 
coefficients before (A — ik)~ 2 and (A — ik)~ 1 leads to the algebraic relations: 

MCM* = (3.20) 
f 1 + M + ^y^j CM f C + M + CM f C + = 0. (3.21) 

As before in order to obtain a nontrivial result M is required to be a degenerate matrix, 
i.e. decomposition ( 13. lip holds true. Then relation ( 13.201) is rewritten as 

(m\C\m)=0. (3.22) 

Relation (13 .2 II) in its turn can be easily reduced to the following linear system for 3- vector 
|n) 

C|n)(m|C\ , . , . 
1+ 1 Q 1 C|m) = ia\n). (3.23) 



2 

by introducing some auxiliary real function a. That linear system allows one to express 
|n) through (m\ and a, namely: 

j^^^-^cj C\m). (3.24) 

In order to find \m) and a we turn back to the equations (I3.15p . Vanishing of the second 
order poles in (13. 15[) leads to the conclusion that 

\m(x, t)) — ipo(x, t, — W)|too), (3.25) 

where |mo) is a constant nonzero 3- vector. After substituting (13.251) into (13. 22[) and taking 
into account ( 12.31) one convinces himself that the components of the polarization vector 
| mo) are no longer independent but satisfy the constraint: 

(m |C|m ) = <^> |m ,i| 2 = |m , 2 | 2 + |^o,3| 2 - (3.26) 

The vanishing condition of the first order poles leads to some differential constraint 
on a(x,t) which is integrated to give: 

o~(x,t) = —K(mo\ip~ 1 (x,t,iK)ip (x,t,iK)C\mo) + o~ , (3.27) 

where <To G M. is a costant of integration. 

Thus to calculate the soliton solution itself one just substitutes the result for \n) and 
\m) into M and uses formula (13. 14ft . As it is seen the new solution is parametrized by the 
polarization vector |mo), the real number <7o and the pole ik. All this can be formulated 
in the following manner: 
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Proposition 2 Let there be given a solution L^(x,t) to \2.13\) , a common solution 
i/jo(x,t,\) to $3. II) . real numbers k > 0, a® and a complex nonzero vector \itiq) satis- 
fying ^3. 2(A) . Then the function L^\x,t) determined by \3.14\) , \3.ll\) . {3.21$ , $3.25\) and 
\3.2Tty is a solution of the system {2.13\) too. The solution ipi(x, t, A) of the dressed linear 
system K3~2}) is given by ipi = ®ip where $ is defined by l[3TQ) . I[J7T1\) . §3~2lj\} , K3~2~3\) 
and ( fOTp . 

One can apply the dressing procedure repeatedly to build a sequence of exact solutions 



-% -% 



*iv T (N) 



More precisely this alternative procedure will be explained in Section 4. 



(3.28) 



3.2 Soliton Solutions 

Let us apply the dressing procedure to the following seed solution 

Lf\x,t) = 




(3.29) 



of equation (12.131) . In this case a fundamental solution to ( 13.1 ft reads: 
ijj (x,t,X) 



( cos(Aa; + fi(A)t)e °s ) 





\isin(Ax + fi(A)t)e _Q ^ — 



isin(Ax + f!(A)t)e^\ 

-2if()(A)t 

e 3 

cos(Ax + fi(A)t)e^ / 



(3.30) 



We recall that fo(A) and fi(A) are even and odd part of the dispersion law induced by the 
Z 2 grading of st(3), see (l2^6|) . 

We are going to consider the generation of a quadruplet soliton first. In this case one 
uses factor (13.81) . It is convenient to decompose the polarization vector |m ) according to 
the eigensubspaces of the endomorphism ip (13.301) : 



\mo) 





+ 7 




(3.31; 



where a,/3,7 are arbitrary complex constants. 

If the vector |mo) is proportional to one of the eigenvectors of the endomorphism ipQ, 
then the corresponding matrix M does not depend on the variables x and t (due to the 
projective nature of the vector |m)) and the corresponding solution (13. 14f) is a simple 
unitary rotation of the constant solution L^. 

Thus elementary solitons correspond to vectors \m ), belonging to essentially two- 
dimensional invariant subspaces of ipo, i.e. they correspond to polarization vectors with 
only one zero coefficient in the expansion (I3.3ip . Let us consider each of these three cases 
in more detail. 



12 



Figure 1: Real and imaginary parts of the component V\ in 13.331 as a function of x. Here 
k=1,u= 1(T 3 , a = 1, p = 1 + i. 



Case (i): a ^ 0, (3^0, 7 = 
The 1-soliton solution is given by: 

u±(x, t) = 

» fx - exp U arctan ^(2^ + 2^ + 0^-0^ 1 (3.32) 
- exp |4i arctan (^ cosh(2ra + 2f , (/ ^ + ^ |a//J|) j j , 

where a = arga, 0, = arg /3. ^(A) and f/(A) are the real and the imaginary part of the 
polynomial fi(A) (resp. f(f (A) and fg(A) stand for the real and imaginary part of fo(A) to 
be used later on). If the dispersion law of NLEE is an even polynomial (fi(A) = 0) then 
the solution f)3.32p becomes stationary: 

ui(x, t) = 

K COs(2o;X + a — 0/3 ) \ 1 (3.33) 



v±(x, t) = exp < 4i arctan 



oj cosh(2/«a; + In \a/j3\) 



A plot of that solution is presented on FigJTJ It is easy to check that u = 0, v = exp(if(x)) 
is an exact solution of ( I2.17P for any differentiable function f(x) tending to when 
x — > ±00. This resembles the case of the three- wave equation [33] where one wave of an 
arbitrary shape is an exact solution of the system and the two other waves are identically 
zero. The solution ( I3.33P has a simple spectral characterisation and an explicitly given 
analytic fundamental solution of the corresponding linear problem. 

If the dispersion law contains odd powers of A as well then the elementary soliton is 
no more stationary. For example in the case of equation (12.181) it reads: 



Ccub(-£; t) 



Ui(x,t) = 0, Vi(x,t) = exp (4i arctan C cub (x,t)) 

Kcos2a;[a; + 8{3k 2 - u 2 )t + (0 Q - <f>p)/2u}\\ (3.34) 



wcosh2/t[x + 8(k 2 - 3u 2 )t + In \(x/P\/2k] 
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Case (ii): a ^ 0, = 0, 7^0 

In this case the solution looks as follows: 



ui(x,t) 



4iuKQ* gen exp i{ux + (f*(ji) + f?(ji))t + <p a - <M 



Sum 



( "- iK)0 - (3.35) 



vi(x,t) = 1-7 . s n2 

where 4> a = arga, 7 = arg7 and 

Q gEn = 2ue KX+{f( ' (p)+ ^ (p))t+lnla/ ~ /l + ( w + i K )e" K:r '" {f o (M)+flI(M)) ^ ln|Q/71 . 

In particular, when f(A) = —A 2 /, i.e. fo(A) = —A 2 and fi(A) = hold, we obtain a 
solution to f l2TTj) : 



Ui{x,t) 



AiumQ* exp i{ujx + (k — uj )t + <p a — </> 7 } 



= 1 



2 (W " iK)g2 (3.36) 



(uj — ik)Q 



where 

Q = 2uje K{x - 2uJt)+lnlahl + (u + i K ) e ~ K ^- 2w *)- ln l a /Tl. 

Contour plots of |wi| 2 and |t>i| 2 of the solutions ( I3.36|) are shown on Figure [2j 

When the dispersion law is odd, say fi_(A) = — 8A 3 , the quadruplet solution represents 
a traveling wave of the form: 

. . 4iu kQ* exp iu\x + 8(3k 2 -u 2 )t+ (d) a - d)^)/uj} 

Ui{X,t) = 



Hojk 

Vi{X,t) = 1 



(to — m)Q 2 

K ^ (3.37) 



(uj — ik)Q 2 
where 

Q = 2u;e K(:r+8(K2 ~ 3tj2) * +ln|a/7l/K) + [uj + j K ) e - K ( a: + 8 ( K2 - 3w2 )*+ ln l Q /7l/«) i 
This is an elementary soliton for the cubic flow NLEE (|2.18p . 



Case (iii): a = 0,/3^0,7^0 

The solution now can be obtained from the solution in the case (ii), by changing a — > 
and x — > —x. 

In the cases (ii) the solution (I3.36P is a soliton of width 1/n moving with velocity 2uj. 
The corresponding soliton in the case (iii) moves with a velocity — 2u. 

In the generic case, when all three constants are non-zero, the solution represents a 
nonlinear deformation of the above described solitons. For k > it may be viewed as 
a decay of unstable time independent soliton from the case (i) into two solitons, corre- 
sponding to the cases (ii) and (iii) (see fig [2}. For k < 0, the solution is a fusion of two 
colliding solitons into a stationary one. 

Let us now consider dressing by a factor with two imaginary poles (doublet case), i.e. 
fi = in. There are two essentially different cases. 
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Figure 2: Contour plot of |tii| 2 (left panel) and \v\\ 2 (right panel) for a generic soliton 
solution (I3.36P as a function of x and t where a = ^ = K = u = l. 



Case (i): a ^ 0, (3^0, 7 = 



From (13.261) it follows that |m ,i| = |m ,3|. It suffices to pick up m ,i 
component is m ,3 = exp(iip), <p G 1R. The doublet solution reads: 



1 and the third 



-gen 



[x,t) 



ui(x,t) = 0, v\(x,t) = exp(4i arctan 5 gen (x,t)) 

o"o — 2k(x + fi (iK.)t) sin ip 
cosh2(/tx + f/(ift)t) + sinh2(/tx + f/(i/t)t) cos<y9 



(3.38) 



If the dispersion law of NLEE is even polynomial, i.e. fi(A) = 0, the 1-soliton solution 
becomes stationary: 



Ui(x, t) = 
v i (x, t) = exp <{ 4i arctan 



ctq — 2Kxsinip 



cosh 2kx + sinh 2kx cos 99 



(3.39) 



Figure 13] presents the argument and the imaginary part of t>i(x) in the stationary case as 
functions of x and the phase </?. 

As in the quadruplet case if the dispersion law is an odd polynomial the doublet solu- 
tion is time-depending. Let us consider the simplest example fi(A) = — 8A 3 corresponding 
to equation (12.181) . Now (I3.38j) obtains the form: 



Ui(x, t) = 



Vi(x,t) = exp(4i arctan E cuh (x,t)) 
o"o — 2k[x + 24K 2 t) sin 99 



(3.40) 



cosh2/t(a; + 8n 2 t) + sinh2K(x + 8n 2 t) coscp 



Case ii. Generic doublet 

Now let us assume mo,2 7^ 0. For simplicity we fix mo,2 = 1- Then the norms of mo,i and 
m ,3 are interrelated through 

I 1 2 I |2 1 

\mo,i\ - m 0)3 = 1- 
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Figure 3: Plots of the argument (left panel) and Imtii(x) (right panel) for the stationary 
solution (13.391 ) function of x and <p; k = o = 1. 



This is why it proves to be convenient to parametrize them as follows: 

m ,i = cosh0 o e i( ^ + ^, m , 3 = I sinh0 o |e i{w -^, (3.41) 
where 0q, <po and (p are arbitrary real numbers. Then the doublet soliton solution reads: 



2A* 

ux(x,t) = __ e i(fo(iK)i+w,) [sinh0+cos^ + isinh0_ sin$ 

, s 2(2io- - 1) 4kr(icx - 1) 

vi(x, t) = l+ A + \ 2 \ 



(3.42) 



where 



A(x,t) = cosh 2 0+ cos 2 (f + cosh 2 0_ sin 2 <£> — ia 

a(x, t) = cr + ftfj (ift)t + k ^a; + fi(W)t^ sinh 20 o sin2y3 

0±(a;,t) = + f 1 / (ift;)t ±0 O - 

Let us consider the special case when the dispersion law is —A 2 /. The solution (I3.42[) 
is significantly simplified if in addition one assumes that m ^/m 01 > {(p = 0). The 
result reads: 

(cosh 2 (ra + ) - i(ab - 2/« 2 t)) 

/ cos1i 2 (kx + O ) + i(cr Q - 2/t 2 t) \ z 2 



(3.43) 

2 (/tx + o ) + i(a -2K 2 t) \ 
cosh 2 (ra + e ) - i(a - 2kH)J (cosh 2 (/ta; + O ) - i(a - 2kH)) a ' 

A plot of Re Ui(x,t) and Re vi(x,t) is shown on Fig. [4] 

It proves to be of some interest to consider the odd dispersion well. In the 

simplest nontrivial situation when fi(A) = — 8A 3 (equation (j2.18j) ) we have 

ui = V V -{ e lw sinh(Ko; + 8kH + O ) 

(cosh (kx + 8k 3 £ + 0n) — icro) 

2 ( 3 - 44 ) 
/ cosh 2 (/tx + 8/? 3 t + O ) + io-o \ 2 



\cosh («a: + 8kH + O ) - i<r / (cosh 2 (/tx + 8/t 3 t + O ) - ia )' 
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Figure 4: Contour plot of Re Ui(x, t) (left panel) and Re vi(x, t) (right panel) for doublet 
soliton ( I3.43P as functions of x and t. Here k = 0, a = 5 and 9 = 0. 

We have assumed above that (p = 0. 



Remark 1 Let us make a few short remarks on the behaviour of doublet soliton ^3.^3 ). 
First of all it is evident that this is not a travelling wave solution. Moreover, as it is seen 
from Fig. [5| the component \u\{x, t)\ 2 has two symmetric maxima and one minimum at 
the origin (resp. \vi(x,t)\ 2 has two symmetric minima and one maximum at the origin). 
The value of the maximum of \ui(x,t)\ 2 (resp. the minimum of \vi(x, t)\ 2 ) first increases 
with time (a if) > 0) and then decreases (o~(t) < 0). The maxima positions of U\ depend 
on t according to: 



(3.45) 



where o~(t) = a — 2K 2 t. The soliton velocity v := d^/dt is not constant but changes 
with t as given by: 



v(t) 



2n 2 ta(t) yi + a 2 (t) 



a 2 (t) 



1 + 



a 2 (t) 



(3.46) 



Such behavior resembles the boomerons and the trappons ^ [5]/. On Fig. it is plotted 
the t- dependence of the soliton velocity. 



3.3 Multisoliton Solutions 

As we have already mentioned the dressing procedure can be applied several times conse- 
quently. Thus after dressing the 1-soliton solution one derives a 2-soliton solution, after 
dressing the 2-soliton solution one obtains a 3-soliton solution and so on. Of course, in 
doing this one is allowed to apply either of dressing factors ( 13. 8 j) and ( 13. 19[) . Therefore the 
multisoliton obtained will be a certain combination of quadruplet and doublet solitons. 
Another way of derivation the multisoliton solution consists in using a dressing factor 
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Figure 5: Contour plot of \ui(x,t)\ 2 (left panel) and \v±(x, t)\ 2 (right panel) for doublet 
soliton solution (I3.43P as a function of x for several values of t: t — 0, 1, 5, 10, 20, 40. 



with a proper number of poles: 



i + aW 

k=l v 



+ 



CM k C 



X- li k X + fi k 



N 2 



i=i 



Pi CPiC 
— f — I f — 

A — iK; A + i/Cj 



(3.47) 



As it follows from ( 13.471) the multisoliton solution obtained will be a mixture of N% quadru- 
plet solitons and N 2 doublet ones. In order to determine the residues of $ one follows 
basically the same steps as in the case of a 2-poles dressing factor. Firstly, the identity 
= II implies that the residues of $ and $ _1 fulfill some algebraic restrictions. For 
example, the condition 



lim (A - fik)®^ 1 = 0, 



k = l,...,Nx (3.48) 

for vanishing of the residue of $<3> _1 at A = jU* leads to the following algebraic restrictions: 

CP/C 



Ni 



r=l 



Ml 



CMjC \ ^ 

7 +^2^ 

r/ 1=1 



fl k + IKi H k - \Ki 



0. (3.49) 



Apart of this type of constraints we have another one originating from vanishing of the 
coefficients before the imaginary poles: 



lim (A - iK,) 2 $$ -1 = {jKifP&P} = 0, I = 1, . . . , N 2) 

A— HKi 

lim d\[{X - iK,) 2 $$ _:l ] = i^©/CP/C + iK,P,C8}C = 0, 

A— s-i/ij 



(3.50) 
(3.51) 



where 



e, = i + k ( Vh — - 



CM fc C 

fi k IKi + fJL k 

| CP S C 

i(Ki + K s 



+ Pl 



CP,C 
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Figure 6: The soliton velocity v(t) and position of the maxima £o(£) of solution (13.431) as 
a functions of t. Here k — a — 1, 9 — and (p — 0. 



Vanishing of the rest of poles of leads to algebraic constraints which coincide with 

(I3.49p - fl3.5ip due to the action of Z 2 reductions. 

Since M k (x,t) and Pi(x,t) must be degenerate matrices one introduces their factor- 
izations M k = \n k ){m k \ and Pi = \qi)(pi\. Substituting it into (I3.49p - p.5ip we reduce the 
first and the third constraint to linear systems for \rik) and \q{) 



Ni 



\m k 



r=l 
Ni 

E 

r=l 



1=1 

N 2 



(3.52) 



where the matrix coefficients read 

(m r \m k ) (m r \C\mk) 
fir ~ fil + fJ-t 

(p s \ m k) (Ps\C\m k ) 



^>rk 


■= /4 ( 


^sk 


■= fi( 




■= -m 


$ si 


■= Ki (- 



c 



IK, 



l»k 



1K S + [I* 



c 



(m r \C\pi) (m r \pi) 
iKi - n k mi + n k 

(Ps\C\pi) (Ps\Pl), 



c 



K s Ki 



K s + Ki 



ft 



S ^ I 



■ (Ps\Ps) ^ 

icr, — C 



By inverting the linear system (I3.52p we can express \n r ) and \q s ) through all \m k ), \pi) 
and &1 and that way determine the dressing factor in terms of the latter. The vectors 
\rrik) and \pi) as well as the functions oi can be found from the natural requirement of 
vanishing of the poles in (I3.15p . The result reads 



\m k (x,t)) = ip {x,t,fi* k )\m kfi ) 
\pi{x,t)) = ipo{x,t, -iKi)\pi fi } 

(Ji(x,t) = -Ki{pi fi \i)~ l (x,t,lKi)i)o(x,t,iKi)C\pi fi ) +<Jlfl. 



(3.53) 
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Analogously to the 2-poles case the components of \pi) are not independent. As a 
result of (13.50!) that the following relations holds true: 

( Pl (x,t)\C\ Pl (x,t)) = (p,,o|Cb,,o) = 0. (3.54) 

Thus we have proved that the dressing factor in the multiple poles case is determined 
if one knows the initial fundamental solution ipo(x,t,\). The multisoliton solution itself 
can be derived through the following formula 

L?(x,t) = <S>(x,t,oo)Lf\x,t)&(x,t,oo), (3.55) 

where 

$(x, t, oo) = l + ^(M fc + CM k C) + ^(P, + CPC). 

k=l 1=1 

From all said above it follows that the algorithm for obtaining the multisoliton solution 
can be presented symbolically as follows 

4 0) — KK)r£i, {|w>}£. ► 

{|^)}£ x , {!»>}£ — ► {M fc }£ 1; {P,}£ — ► L« 



4 Interactions of Quadruplet Solitons 

In this section we aim to study the interactions of solitons we have derived. We shall 
restrict ourselves with quadruplet solitons for NLEEs with odd dispersion laws. This is 
the simplest case since the solitons are travelling wave-type solutions. The interactions of 
the other types of solitons require a special treatment and will be done elsewhere. 

Our study will be based on the Zakharov-Shabat scheme [31] applied to the recur- 
sive procedure (I3.28p . Their approach consists in calculating the asymptotics of generic 
7V-soliton solution for t — > ±oo and establishing the pure elastic character of the inter- 
actions of generic soliton, i.e. solitons travelling at different velocities. The pure elastic 
character of the soliton interactions is demonstrated by the fact that for t — > ±oo the 
iV-soliton solution splits into a sum of iV one soliton solutions preserving its amplitudes 
and velocities. The only effect of the interaction consists in shifting the center of mass 
and the initial phase of the solitons. 

The 1-soliton dressing factor corresponding to the quadruplet case with poles at ±fi k 
is given by: 

$ k (x, t, A, Hk) = 1 + — ^M fc (x, t) + -^-CM,(i, t)C. (4.1) 
A - fik A + fi k 

The residues M k (x,t) = \nk)(m k \ are determined by the following equalities 

K) 

\m k (x,t)) 

Let us now outline the alternative procedure for constructing the iV-soliton solutions of 
the NLEE (I2.13p . The idea is to apply subsequently N times the the one-soliton dressing. 



1 / (m k \m k ) (m k \C\m k ) 



ip (x,t,fi k )\m k0 ), 




(4.2) 
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For simplicity we assume that all N solitons are of quadruplet type. As a result the 
sequence of mappings ( 13. 28ft allows us to constructs a sequence of Lax operators with 

(k) 

potentials L\ , k = 1, . . . , N and eigenfunctions: 



(4.3) 



xf k) (x,t,\) = & k (x,t, A,/xjfc)$ fe _i(x,t, A,/x fe _i) . . .$i(a;,t, A,/ii) 

x VoOM, A)$ t lj _(A, / u 1 ) . . . ^_ 1 _(A,/i fc _i)^ _(A,/i fe ), 

where 

$ fc -(A,/ifc)= lim * fc (a?,t,A,^ fc ) (4.4) 
The dressing factors $> k (x,t, A,^) are constructed in analogy with ( 14. ip as follows: 



(4.5) 



$ fc (x, t, A, (i k ) = TL + — ^M fc (rx, t) + —^CM k (x, t)C 

A - jJLk A + yUfc 

M k (x,t) = — V 1 7 - 1 1 7 C m fc (mJ 
1 ' J nl\ 2in k 2u k ) 1 fe| 

|m fc ) = $fc-i(x,t, A*fe=i) • • • $i(x,t,/i fc ,/ii)|m fe ). 
Thus for the iV-soliton potential we obtain: 

L[ m (x,t) = limxf N) {x,t,\)L ( f } .xt N) {x,t,\) (4.6) 

Next we recall that we are considering NLEE with odd dispersion laws ( 12.13b ). 
Their one-soliton solutions are traveling waves and depend on Z k = x — V k t, where 
V k = l/KfcIm/i(/i fc ). In particular, for the eq. ( 12.181) fi(A) = — 8A 3 and V k = 8(3/4 — k\). 
Now let us to pick up the trajectory of the iV-th soliton: Z N = x — 2u)Nt/3 = fixed and 
evaluate the asymptotics of L^(x,t) for t — > ±oo for fixed Z N . This will allow us to see 
what are the effects of the soliton interactions on the N-th soliton. 

In what follows we will assume that all solitons move with different velocities, i.e. 
Vj 7^ Vfc for k 7^ j. It is natural to split the solitons in two groups: 

M + = {V k : V k > V N }, M_ ee {V k : V 3 < u N }, (4.7) 

i.e., the solitons belonging to M + are moving faster than the N-th soliton, while the ones 
belonging to M_ are slower. 

Now we are able to calculate the limits of $fe(x, t, A) for t — > ±oo for fixed Z^. To do 
this we firstly need to obtain the limits of the one-soliton dressing factor for x — > ±oo. It 
can be verified that: 

Cfc(A) 0-c' fc (A), 
$fc,-(A,/i fc ) := lim $ k (x,t,\)=\ 10 

-4(A) c k (X) 

■ Cfc (A) 4(A) 
$ fej+ (A,/i fc ) := lim $ k (x,t,\) = I 10 



(4.8) 



where 



4(A) Cfc (A) 



m _ A 2 - \Hk\ 2 f^k Kf^k — f4s) 
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Note that the asymptotics $fe,±(A, fi k ) do not depend upon the polarization vectors \rriko) 
and that they commute for different values of A. This allows us to describe explicitly the 
iV-soliton interactions of quadruplet solitons. 

The action of § kt ±(\,fi k ) on the polarization vectors produces the equalities: 

fa k + P k \ /af + ^f 
$k,±(A,/i*) Ik = Ik 

\a k -(3 k J \a£-p£/ (4.9) 

= fj, k X±fj,% £f = n k A =f /4 
a k yL* k \±yi, k (3 k /ifcAT/^fc' 

Next we have to evaluate the asymptotics of \m k (x,t)) when t — > ±00 along the 
trajectory Z^{x,t) = const. This is done recursively using (14. 8p . Skipping all technical 
details here we get: 



\m N (x,t)) ~ IT $+(/i iV ,// J ) I f ^-(fi N ,fXj)\m N (x,t)) 



t— >oo 

jeM+ jeM- 



\m N (x,t)) ~ IT ®-(n N ,Hj) TT ^ + (fjL N ,fjLj)\m N (x,t)). 



(4.10) 



t— 00 

jeM + jeM- 



Then from (14. 9p and (I4.10p one deduces that: 

+ N N 

11 . 11 11 - : ir. 11 11 > 



N,j 

fc=i ™ jeM_ fe=i ^ k j&A+ jeM- 

fl+ N Ft- N 

Pn TT ^ TT d TT 1 P. 



; 11.. 11 11 ■ f=n§n^n^ (4ll) 

UN + /i* fJ-N — f*j 

A N,j — j , N,j — • 

/ijV + fJ'N ~ A*j 

As a result we obtain that: i) the soliton interactions are purely elastic, and ii) their effect 
is shifts of the relative center of mass and the phase 5n = arga — arg/3 of the solitons: 

5^ = 5jv ± ^ 0^ T ^ ( 4 - 12 ) 
%N)j = 77— ( m l^ivjl - In l-Bjvjl), <?W = arg(v4 A r J ) - arg(S iV)i ). 

5 Integrals of Motion 

Here we will sketch briefly the direct method for finding integrals of motion, introduced 
by Drinfel'd and Sokolov [6j. We will apply it to the system (I2.17p . In order to do that it 
proves to be technically more convenient to deal with the Lax pair (12. 19f) , (12.2ip . We will 
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use the transformation 7(x, t, A) that diagonalizes simultaneously the Lax pair L and A: 

f\ 

L = y- l L7 = id x + A J + L + -i + • • • 

A-i (5 ' 1} 
A = ?- l A? = \d t + A 2 / + AA_ a + A + - 1 + • • • . 

A 

Here all matrix coefficients -A-i and Ak, k = 0, 1, . . . are diagonal. Using the asymp- 
totic expansion for 7(x,t, A): 

y(M,A) = * + ^ + ^ + --- (5.2) 
A A 



one can get a set of recurrence relations: 



Uo + Jpi = Lo+PiJ (5.3) 
ipi,x + U pi + Jp 2 = Li+ Pi£ + ViJ (5.4) 



fc-i 



ipk,x + U p k + Jpk+i = £fc + Pk+i J + 2^ 

Pk—m^m (5-5) 



m=0 



Here we assume that all coefficients pi (1 = 1,2,...) are off-diagonal matrices. 

In order to solve the recursion relations above, we will split each relation into a diagonal 
and off-diagonal part. For example, treating this way the first relation above one gets 

L = U$, Ul = -[J, Pl ], (5.6) 

where the superscripts d and / above denote projection onto diagonal and off-diagonal 
part of a matrix respectively. Taking into account the explicit form of Uq for £q we have 



1 
0-2 
1 



L = X -(uul + vvl) \ 0-2 J . (5.7) 



Thus as a density of our first integral we can choose: Jo — u*u x + v*v x . It represents 
momentum density of our system. For the stationary solutions (13.331) and ( 13.391 the 
momentum density is depicted on Figure [71 It is evidential, that the integrals of motion 
are well localized function of x. 

Similarly, for the second integral density one gets 



Ji = \uu* x + vv* x \ + A\uv x — vu x \ . 
In general, the-k integral of motion can be calculated through the formula 

L k =(yi Pk y. (5 A 

The matrix p^ in its turn is obtained from the following recursive formula 



fc-i 



p k = -ad j 1 \Pk~l,x + (U Pk-l) f - Pk-l-m^m ■ (5.9) 



m=0 
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Figure 7: Plots of the density of the first integral of motion as a function of x evaluated 
on the stationary quadruplet soliton ( I3.33|) fora = /3 = 7 = w = l,<5 = (left panel) and 
for the stationary doublet soliton (I3.39j) for k = a = 1, = (right panel). 

Note, that the zero curvature representation is gauge invariant, i. e. [£,.A] = is 
fulfilled. Since [£fc,./bj = the commutativity of £ and A is equivalent to the following 
requirements 

d x A- 1 = 0, d t L k -d x A k = 0, A; = 0,1,... (5.10) 
Hence represent densities of the integrals of motion we are interested in. 

6 Conclusions 

The soliton solutions for a hierarchy of NLEEs related to the symmetric space SU(3) / S(U(1) 
xU(2)) are constructed. In order to obtain the soliton solutions we have applied the dress- 
ing procedure with a 2-poles dressing factor. It has been shown that there exist two types 
of 1-soliton solutions: quadruplet solitons which are associated with 4 symmetrically lo- 
cated eigenvalues of L and doublet solitons which are associated with a pair of purely 
imaginary eigenvalues. This remarkable fact is a consequence of the simultaneous action 
of two Z 2 reductions on the Lax pair. The properties of the elementary solitons depend 
crucially upon the symmetry properties the dispersion law. For example, if the dispersion 
law is an even polynomial then the elementary soliton of the first type will be stationary 
(see formula (13.331) ) otherwise it is time-dependent (formula (I3.34p ). In the case of the 
doublet type solitons the situation changes significantly — the components of the polar- 
ization vector \m ) are no longer independent, see (13.261) . This is why we have only two 
cases possible: generic case and a degenerate case. In the latter case the doublet soliton 
is stationary if f(A) is an even polynomial, otherwise they are time-depending. In the 
generic case a new phenomenon arises. When the dispersion law is an even polynomial 
the soliton is not a traveling wave. Its behavior resembles that of trappons and boomerons 
— the soliton velocity is not fixed but varies with time. 

We have described the quadruplet soliton interactions for NLEE with odd dispersion 
laws by calculating explicitly their asymptotics along the soliton trajectories in the generic 
case (different soliton velocities). The important result consisted in the following: 

i) the iV-soliton interactions are purely elastic and always split into sequences of ele- 
mentary 2-soliton interactions; 

ii) the effect of each 2-soliton interaction consists in shifts of the relative center of mass 
and relative phases of each of the solitons; 
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iii) the corresponding shifts are different from the ones for the NLS and Heisenberg 
ferromagnetic equations. 
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